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We simulate a growth model with restricted surface relaxation process in d = 1 and d = 2, where 
d is the dimensionality of a flat substrate. In this model, each particle can relax on the surface to a 
local minimum, as the Edwards- Wilkinson linear model, but only within a distance s. If the local 
minimum is out from this distance, the particle evaporates through a refuse mechanism similar to 
the Kim-Kosterlitz nonlinear model. In d = 1, the growth exponent j3, measured from the temporal 
behavior of roughness, indicates that in the coarse-grained limit, the linear term of the Kardar-Parisi- 
Zhang equation dominates in short times (low-roughness) and, in asymptotic times, the nonlinear 
term prevails. The crossover between linear and nonlinear behaviors occurs in a characteristic time 
tc which only depends on the magnitude of the parameter s, related to the nonlinear term. In d = 2, 
we find indications of a similar crossover, that is, logarithmic temporal behavior of roughness in 
short times and power law behavior in asymptotic times. 
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I. INTRODUCTION 

Self-affine interfaces generated by noncquilibrium sur- 
face growth have been intensively studied in recent yeacs 
[1-3] . Kinetig roughening models as ballistic depositionQ, 
Eden modftfl, soHd-on-solid model (SOS) with surface 
relaxationHfl, SOS with refuseH and SOS with diffusionEI 
are some examples of growth models that belongs to dis- 
tinct universality classes. Theoretically, these growth 
processes are studied in three different schemes: through 
a continuum description using Langevin-like equation 
and renormalization group analisys for solving it; through 
numerical solutions of these equations; and through com- 
puter simulations of discrete models. The main goal is 
to obtain the universality class of a specific model and to 
get information about the presence of nonlinearities and 
broken symmetries. 

In computer simulation of lattice growth models, inter- 
faces are described by a discrete set {hi(t)} which repre- 
sents the height of a site i at the time t. Such an interface 
has L'^ sites, where L is the linear size and d is the di- 
mension of the substrate. The roughness of the interface 
uj is defined as the root mean square of the {hi — h} 
distribution, 

c^HL,t) = /^i-^j2(^,-hry (1) 

where h is the mean height at time t and the angular 
brackets means the average over independent samples. 

The universality class of a discrete growth model is 
obtained through the temporal and spatial behaviors of 
roughness. In the most of kinetic roughening processes, 
which starts at t = from a flat substrate, the temporal 
behavior of roughness is described by the power law be- 
havior, uj{L,t) ^ , when 1 <C i <C tx , and its spatial be- 
havior, in the steady state, is described by ujsat{L) ~ L", 
for t^tx- The exponents P and a are the growth and 
roughness exponents, respectivelly, and tx is the satura- 
tion time. These two behaviors are joined at the Family 
and Vicsek dynamical scaling relationtB 

u;iL,t) ^ Ly (2) 

where the function f{x) must be L-independent. This 
function scales as f{x) ^ for short times and tends to 
a constant in the steady state. The dynamical exponent z 
is related with a and /3 through the relationship z — a/ (3 
and it shows how the saturation time depends on the 
system size L: tx ~ L^- Two of these exponents predict 
which universality class a model belongs to. 

These universality classes are related to the dominant 
term of the stochastic differential equation in the contin- 
uum limit. In the next section, we describe two stochastic 



differential equations which represent two distinct univer- 
sality classes. The first is a linear equation proposed, in 
1982, by Edwards and Wilkinson (EW equation)^. The 
second is a nonlinear equation, introduced, by Kardar, 
Parisi and Zhang in 1986 (KPZ equation)!!!. When the 
nonlinear term of this equation is null, the EW linear 
equation is recuperated. 

In this article, we report on simulations of a growth 
model with restricted relaxation process (called RR 
model), which was proposed in order to study the 
crossover between the linear and nonlinear regimes of 
the KPZ equation. The RR cotmbines features of the 
Edwards- Wilkinson (EW) modeHfl, related to the EW 
linear equation, and the Kim-Kosterlitz (KK) modelB, 
related to the KPZ nonlinear equation. The crossover be- 
tween linear and nonlinear regimes of the KPZ equation 
was studied numerically in d = 1 and d = 2 through nu- 
merical solutions of this equation where variations of the. 
amount of nonlinearity was allowed. Moser and KerteszEj 
found, in d = 1, the growth exponent (3 = 1/3 for 
all values of the nonlinearity. In d = 2, the authors 
have fouftd (3 = 0.240, close to the Kim and Koster- 
litz valuecl. So, a crossover between linear and nonlinear, 
regimes was not verified. Grossmann, Guo and GrantEj 
also have obtained numerical solutions of KPZ equation 
where the surface tension was fixed and the nonlinear 
parameter was continuously changed. For small systems, 
the authors have found growth exponents in the inter- 
val 1/4 < (3 < 1/3. So, the crossover was character- 
ized by a continuum change in the growth exponent, 
which in c? = 1 means that the dynamical exponent z 
depends on the amount of the nonlinearity. This contin- 
uous crossover |iH |6L= 1 was also verified by simulations on 
growth mpdelgliO. Theoretically, in d = 1, Nattermann 
and Tang|!3 have studied the KPZ equation in the low 
nonlinear limit through renormalization group analisys 
and they have obtained a different result: two behaviors 
separated by a characteristic time tc- For t <^ tc, the 
linear behavior of roughness was found, while for t ^ tc, 
the nonlinearity dominates. 

This article is organized in sections. The next section 
presents the general approach for kinetic roughening and 
describes the EW and KK models. In section III, we in- 
troduce the RR model and we show our results which cop-, 
firm the theoretical previsions of Natterman and Tangl!3 
and, in section IV, we finally show our conclusions. 

II. THEORY AND DISCRETE MODELS 

In the continuum (coarse-grained) description, the in- 
terface motion is described through Langevin-like equa- 
tions [1-3], 

^^^=^;o + r/(x,^) + $[Mx,^)] . (3) 
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In this equation, vq and 77(x, t) are the deposition rate 
and its noise, respectively. This white noise has zero 
mean and variance given by 



(r]{yi, i)7/(x',t')^ = D6'^{yi - x')6{t - t). 



(4) 



$ [/i(x, t)] is the term related to the correlations between 
neighbors which can have linear and nonlinear functions 
of /i(x, t). 

Our intagst is focused on the SOS model with surface 
relaxationQ'EI and at the SOS model with restriction^. 
The SOS model with surface relaxation, introduced by 
Edwards and WilkinsonO in 1982, is a random deposi- 
tion of particles where the difference of height constraint 
between the neighbors {j} of a site i is given by 



ho) < M, 



(5) 



where Al is the parameter that controls the roughness. 
In this work, we always use M — \. If the height of 
the deposited particle on the site i does not satisfy the 
height constraint, this, particle must be moved to a lo- 
cal miminum. FamilyQ have obtained for this model, in 
d = 1, the exponents /3 = 0.25(1) and a = 0.48(2). This 
result indicates that this model, in a coarse-grained limit, 
belongs to the universality class defined by the EW equa- 
tion 



dh{x,t) 
dt 



= wo + ??(x,i) + i/V2/i(x,i), 



(6) 



where the laplacian term is related to the surface relax- 
ation process. The exponents ofrthis linear equation, ob- 
tained through Fourier analysisBO, are /3=l/4, a=l/2 
and z = 2 for d = 1. For d = 2, which is the critical di- 
mension of Eq.6, these exponents are P = a — which 
means that the roughness has logarithimic behavior in 
space and time (tj^ ~ logt, for t ^ L^, and w^^jj ~ logi, 
for t 3> i^). Both, the EW model and the EW equation, 
generate gaussian height distributions as well. 

In the SOS model with restriction, particles are also 
randomly deposited onto a substrate and the difference 
of height constraint is the same of the EW model, Eq.5, 
but any kind of relaxation is allowed. If the height of a 
deposited particle does not satisfy Eq.5, this choice must 
be refused, that is, the particle evaporates. |-|This model 
was proposed in 1989 by Kim and KosterlitzEl in order to 
study nonlinear kinetic roughening in high dimensions. 
They numerically showed, in d = 1, that this model, 
named KK model, belongs to the universality class of the 
well known— KPZ equation, proposed by Kardar, Parisi 
and ZhangO 

^^^^ = i;o + r7(x, t) + uV^hi^, t) + ^ [V/i(x, t)f . 

(7) 

The appearance of the nonlinear term [V/i(x, t)]^ is due 
to the lateral growth, that is, the dependence of the 



growth velocity on a local normal of the growing inter- 
face, or to the appearance of a perpendicular driven force 
that leads to a growth velocity greater or smaller than 
the deposition rate vq. In the case of the KK model, 
for example, the refuse mechanism makes the growth ve- 
locity smaller than the deposition rate. In d = 1, the 
exponents of this equationU are (3 — 1/3, a — 1/2 and 
z = 3/2. In d = 2, the analytical solution is not kiiown. 
In c? = 1, numerical simulations of the KK modela indi- 
cate P — 0.332(5) and the a-exponent close to the ex- 
pected value (a — 1/2) and, in d = 2, /3 = 0.250(5) and 
a = 0.40(1). 

Eq.7 is not invariant under the h ^ —h transforma- 
tion, which means that the up-down symmetry is broken 
in surfaces generated by a KPZ process. This fact leads 
to deviations in the gaussian caracter of the height dis- 
tributions which can be measured using other moments 
of the distribution. The Eq.l can be generalized for any 
moments of height distribution as 



/ L 



(8) 



where q is the order of the moment. Note that the rough- 
ness uj{L,t) is related to the second moment: uj'^{L,t) = 
W2{L, t). A growing profile has up-down symmetry when 
positive and negative local curvatures are equals, and, in 
this case, W3{L,t) vanishes. On the other hand, when 
asymmetries are present, W3(L,t) ^ 0. The skewness, 
defined by 



SiL,t) = 



W3iL,t) 



(9) 



is the ideal measurement of deviations from gaussian be- 
havior. In the case of the EW model, S{L, t) — Ojalways. 
For systems in KPZ class in d = 1, Krug et aZ.EZI have 
indicated \S{L,t)\ « 0.28 as an universal value in the 
transient state. In the steady state, the profile shows a 
random-walk character, that is, the height distribution is 
gaussian and Ssat{L,t ^ ix) = 0. In d = 2, numerical 
simulations of the KK model indicate S w —0.40, ip-jljfef 
transient state, and Ssat ~ 0.28, in the steady statetSlEJ. 



III. MODEL DESCRIPTION AND RESULTS 

In this article, we report on simulations of the growth 
model with restricted surface relaxation model (RR 
model). In the EW model, each incoming particle 
must search the local minima when the height constraint 
(Eq.5) is not satisfied. We introduce a parameter s that 
is the number of lattice units allowed for the relaxation 
process. If the deposited particle does not find the mini- 
mum after s relaxations, then this choice must be refused, 
as in KK model. For s = 0, the KK model is recuperated 
and s ^ oo yields the EW model. 
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A. d = 1 Results 

Figure 1 shows the log- log plot of roughness uj{L,t) 
vs time t for a system with L = 10^ sites and s = 2. 
The time unity means L'^ attempts of deposition. The 
two straight lines in the figure are showing the power-law 
fits with (3 = 0.249(1) {1 < t < 10^) and /? = 0.332(1) 
(10"^ < t < 10^). The intersection of these two lines 
define the crossover time tc- This crossover is easy to un- 
derstand considering the increase of roughness: in short 
times, when the roughness is still small, particles do not 
need to relax very much for scaching local minima, so the 
linear behavior might dominate. In large times, on the 
other hand, we note the appearance of relaxation lenghts 
bigger than those observed in short times. As relaxation 
processes are linked to refuse processes in this model, the 
system undergoes a crossover to the nonlinear behavior. 
The crossover time tc is independent of the system size 
L and it is only a function of the parameter s. 

For a better understanding of this crossover wc study 
the statistics of relaxations in the EW model, where 
relaxations of all sizes can occur. At the time t, be 
{Nk{t)) the mean number of particles which diffused k 
sites searching local minima. So, k relaxations occur with 
probability 



Pk{t) = 



(10) 



These probabilities have an initial temporal dependence 

and a steady state Pfc(oo), whose values arc shown in 
Table 1. We note a strong decrease of the probability 
with the number of sites diffused, which shows the rare 
ocurrences of the relaxation with large relaxation lenghts. 

It is interesting to analyse the approach of this proba- 
bility to the steady state because each Pk (t) has (lifFc;rcnt 
convergence times. For doing this, we define a normalised 
probability of k relaxations as 



Pk{t) = 



Pkjt) 

Pfe(oo) 



(11) 



Figure 2 shows the plots of Pk{t) vs t for fc = 0; 1; 2; 3; 4 
for the EW model with L = 10^. Note the differences 
among the convergences to each steady state: the tem- 
poral behaviors of po and pi quickly go to its steady val- 
ues, while p2 and p^ tend to unity only at i w 10^ and 
t ~ 10^, respectively. This fact suggests that large re- 
laxation lenghts might occur at large deposition times 
with small probabilities. So, drawing attention to the 
curve P2{t) vs t, we observe that the saturation occurs 
at f « 10"^. In the RR model with s = 2, we estimate 
tc « 10^, which indicates that the crossover from lin- 
ear to nonlinear behaviors of the RR model occurs, for a 
value of s, when pk=s(t) is time independent. 

This behavior of pk is responsible for the dependence 
of the crossover time tc with the parameter s in the RR 
model. Figure 3 shows clearly this dependence with the 



plots oiui{L,t)/t'^ vs t for: (a) (3=1/4: and (b) /? = 1/3, 
with s = 0, 2, 4 and L = 10'"'. In (a), the curve with s = 
always grows, while the curve for s = 4 remains constant, 
indicating that f3 = 1/4, for this value of L, is the correct 
value for s = 4 (EW behavior) and a noncorrect for s = 
(KK nonlinear behavior). In (b) similar conclusions are 
obtained with /3 = 1/3. For s = 2, the initial linear 
and the asymptotic KPZ behaviors are well observed in 
(a) and (b), respectively. Note that the crossover for 
s = 2 occurs when the P2{t) is time independent, that is, 
when P2{t) « 1, in Figure 2. For s = 4, we observe only 
the linear behavior because the total deposition time is 
smaller than the crossover time tc- 

In order to do a more complete characterization of this 
crossover, we also analyse the temporal behavior of the 
skewness S{L,t). Figure 4 shows plots of S{L,t) vs t 
for L = 10^ and several values of the parameter s. For 
s = 0, the skewness S{t) goes quickly to the KPZ tran- 
sient value S = —0.28. As we have explained, S ^ 
means that the interfaces have not up-down symmetry. 
We find for s = 2, a slower approach to the KPZ value, 
than observed for s = 0, indicating that the up-down 
symmetry is gradually lost when s > 0. In particular, for 
s = 3, we clearly observe an initial behavior where 5 ~ 
and an approach faster to the KPZ value at t « 10^. The 
temporal dependence oi S{L,t), which is an independent 
measurement of the universality class, also indicates an 
initial linear behavior and the KPZ behavior in asymp- 
totic times. 



B. d = 2 Results 

It is also interesting to study this crossover in d = 
2, because changes in the morphology are expected in 
this dimension. For the KK model, we have power-law 



0.78^ 



divergencies of roughness (cj ^ f and L' 
while the EW model shows logarithmic divergences {lv"^ ~ 
Int and ~ InL. 

In order to avoid saturation effects, wc work with 
L = 2000 (4x10^' sites), and we do simulations until 
t = 10*. Due to the computational cost, we perform 
only two samples for each value of the parameter s and, 
consequently, the data quality in this section is poorer 
than in the previous subsection. The crossover is anal- 
ysed through a similar process done in d = 1. We expect 
logarithmic behavior when the linear term dominates, so 
we define 



io^{t)-B, 
As ln{t) 



(12) 



where Ag and Bg are the s-dependent coefficients ob- 
tained through a logarithmic regression in the interval 
2 < t < 100. In Figure 5, we plot the temporal behav- 
ior of the function Ys{t), in a semi-logarithmic scale. If 
the temporal behavior of the roughness is logarithmic, 
Ys{t) must be a constant, equal to one. We note that 
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this occurs for s = oo, but, for s = 1, this constant be- 
havior remains until t « 10^. After this time, there is a 
crossover to the power law behavior. However, until the 
studied time, it is not possible to determine the exponent. 

In order to verify the existence of a power law behavior 
with the KPZ exponent /3, we need to analyse values of 
the parameter s smaller than s — 1, that is, we need to 
do continuum variations in s. So, we do this by assinging 
a probability s for the particle relax one lattice unit and 
a probability (1 — s) for the particle be evaporated. In 
figure 6, we show the graph of w/t^ vs t, for s = (KK 
model) and s = 0.1,0.3,0.5. In this graph, it was used 
f3 = A/4 which is the exponent for the KK model in 
d = 23. If the temporal behavior of the roughness has a 
power law behavior with this exponent, the curves must 
be horizontal. For the RR model (s 7^ 0), the curves 
show asymptotic approaches to this nonlinear behavior. 
For greater values of the parameter s, the crossover time 
to this behavior occurs for times greater than the time 
studied (^=10^). 

IV. CONCLUSIONS 

We have studied the model with restricted surface re- 
laxation which combines the main features of the model 
with surface relaxation (EW model) and the model with 
refuse (KK model). A power-law temporal behavior of 
roughness with two growth exponents was observed in 
d — I: The linear growth exponent, /3 = 1/4, occurs in 
short times and the nonlinear one, j3 = 1/3, appears in 
the asymptotic limit. This result suggests the follow- 
ing description: the linear term of the KPZ equation 
(Eq.7) dominates in short times and the nonlinear term 
dominates in asymptotic times. We also noted that the 
crossover time tc is independent of the system size L and 
it is only function of the parameter s. This result cor- 
roborates the renormalization group solution made by 
Nattermann and Tang|l3, where the KPZ equation with 
small nonlinear term was considered. In d = 2, we have 
found indications of the same kind of crossover: A log- 
arithmic temporal behavior of roughness in short times, 
which is related to the linear EW equation, and a power- 
law behavior with /3 — 0.25 in asymptotic times, related 
to the nonlinear KPZ equation. 
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Figure captions 



TABLE 1 



Figure 1 

The temporal behavior of the roughness uj{L, t) for the 
RR model simulated at a substrate with L = 10'^ sites 
and s = 2, in a log- log plot. The two straight lines are 
showing the power-law fit results with (3 = 1/4 (short 
times) and /3 = 1/3 (asymptotic times). The crossover 
time tc is defined as the intersection of these two lines, 
as indicated in the figure. 

Figure 2 

The temporal behavior of normalised probability of k 
relaxations Pk{t) for the EW model with L = 10^ sites 
for A; = 0; 1; 2; 3; 4, from top to below. 

Figure 3 

The log-log plots of temporal behaviors of: (a) uj/t^/'^ 
and {h)w/t^/^ for s = 0: 2; 4 and L = 10^ We might ob- 
serve, for s = 2, when the system drives away from the 
linear behavior in (a) and, the arrival at the nonlinear 
KPZ behavior in (b). For s = 4, the system must ap- 
proach to the KPZ behavior at deposition time greater 
than i = 10^. 

Figure 4 

The temporal behavior of the skewness S{L,t) for 
L = 10^ for several values of the parameter s. The two 
horizontal lines show the value 5 = (EW value) and 
S = -0.28 (KPZ value). 

Figure 5 

The temporal behavior, in a semi-logarithmic scale, of 
the function Ys{t) for s = 1 and s = oo (EW) with L = 
2000. The dashed horizontal line indicates the constant 
value Ys{t) = 1. 

Figure 6 

The semi-logarithmic graph of uj{L,t)/t^^^ in function 
of the time t for L = 2000 and some values of the param- 
eter s. 



k Pfc(^) 

0.5813(3) 

1 0.3481(3) 

2 0.0616(1) 

3 0.00812(2) 

4 0.00085(1) 

5 7(5) X 10-5 



The steady state values of probability of difl[usion k sites 
determined as the mean value in the interval 10* < t < 
10^ for a EW system with L = 10^. 



6 



(a) °s = 
+ s = 2 
* s=4 



fH \ — h 



(b) 



3 Oo o o o ooooooooooooooooooooooooooooo 



t-+^ 



-7|f sL 



10 10 



10 



10 10 



T. J. da Silva and J. G. Moreira - Figure 3 



< > 



0.00 



o 



-0.10 



-0.20 



□ 

* □ 



□ 



□ 

* □ 

* °n 
* □ 



□ r 



* s = 

□ s = 0.1 

• s = 0.2 

o s = 0.3 
+ s = 0.4 



-0.30 



10' 



1& 



10' 

t 



10^ 



10^ 



T. J. da Silva and J. G. Moreira - Figure 4 



1.2 



* 
* 



1.1 



^ * ** 
* *** 

* * 



c 



1.0 



* * 



,oo. 



° ° o Co o 



cP 



O 



10^ 



10^ 



10 



10^ 



T. J. da Silva and J. G. Moreira - Figure 5 



oo 



ooo 



"oooooooooooooooa 



'ooco'- 



X ****** 

******* 



— 



*s = 0.1 'x^. 
* s = 03 



X 

X 



5 = 0.5 



Xx 

Xx 

Xy 

Xxv 



^XXxx 



"XXxX. ^ 

XXXxXxxX 

•^xxl 



' ' 

lO' lO' 10^ lO"" 



T. J. da Silva and J. G. Moreira - Figure 6 



